Landau hydrodynamical model at RHIC and LHC 
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The rapidity distribution and transverse spectra of most copious particles such as pions, Kaons 
and antiprotons from central Au+Au collisions at -^sjvjv = 200 GeV and central Pb+Pb collisions 
at y/SNN ~ 2.76 TeV have been investigated in the framework of Landau hydrodynamical model. 
With a more realistic choice of freeze-out condition and the employment of lattice equation of state, 
we find transverse expansion of the collision systems is important to explain the observed data. 
With the increase of collision energy from RHIC to LHC, transverse flow becomes more and more 
important for hadron production at midrapidity, especially for more massive particle. 

PACS numbers: 



I. INTRODUCTION 

The Relativistic Heavy Ion collider (RHIC) and the 
Large Hadronic Collider (LHC), the biggest high energy 
nuclear experiments in the world, realized the creation of 
a quark gluon plasma (QGP), a thermalized system con- 
sist of free quarks and gluons, virhich was the main matter 
of our early universe shortly after the Big Bang. In order 
to understand the organ of our universe and to study the 
non-perturbative region of quantum chromodynamics, it 
is of great importance to measure the macroscopic prop- 
erties of this matter, such as equation of state, transport 
coefficient, viscosity, phase boundary, and so on. How- 
ever, direct measurement is extremely difficult because 
the QGP matter made-in-lab exists only shortly in vio- 
lent systems, due to the confinement of quarks. In fact, 
those properties can only be extracted from experimental 
data based phenomenological models. 

Among those models, statistical thermal model was 
ffi'st proposed in 1950's pj, though there was neither the 
terminology of quark nor that of QGP at that time. This 
model is still in use today. The basic assumption is, after 
the creation of a fireball (of QGP), all hadron species are 
strongly produced under the same temperature followed 
with weak decays and hadronic interactions. The phase 
transition temperature Tc is one of the main parameters 
of the model, which can be extracted from experimen- 
tal data of particle yields/ratios. The obtained temper- 
ature is about 160 MeV, very close to Lattice (LQCD) 
prediction^]. 

However, statistical thermal model can neither de- 
scribe the dynamics evolution of the collision systems 
nor cover more observables than particle yields/ratios. 
Soon Laudau Q proposed hydrodynamical model, which 
made one step forward. With latter development, more 
hydrodynamical models have been employed in heavy ion 
physics, from Bjorken's one-dimensional boost-invariant 
hydrodynamics to hydrodynamics in two-dimensional 



and three-dimensional space, additonally with viscosity 
and fluctuated initialization. Thus hydrodynamics has 
become one of the most powerful models, with more and 
more experimental observables explained, with a more 
and more realistic but also complicate hydrodynami- 
cal evolution. In this paper we will develop the recent 
work0,Q, investigate the employment of the olddest and 
simplest hydrodynamical model. Landau hydrodynami- 
cal model, to RHIC and LHC experiments. 

The paper is organized as following. In section II, 
we will introduce our calculation approach of particle 
production based on Landau model. In section HI we 
will present our results on both rapidity distribution and 
transverse spectra of the most copious particles such 
as charged pions, Kaons and antiproton from central 
Au+Au collisions at y^s^vTV ~ 200 GeV and central 
Pb+Pb collisions at y^ijvjv = 2.76 TeV, and compare 
with RHIC and LHC data. The effects from different 
freeze-out conditions and equation of states will be shown 
as well. In section IV we will draw the conclusions on 
what we have learned from this exercise. 



II. CALCULATION APPROACH 

Hydrodynamics is often employed to describe the evo- 
lution of the collision system, from an initial time (ini- 
tial condition) to a freeze-out time (FO condition) . Dur- 
ing this period, the system is governed with the energy- 
momentum conservation law 



0. 



(1) 



Solving this fluid equation with a given dynamical 
equation of state (EOS), and decomposing the energy- 
momentum tensor 



T"'' = {e + P)u''u'' - Pgf' 



(2) 
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one can obtain the energy density e, pressure P, flow ve- 
locity u"^ at each space-time point of the collision system. 
Here we assume an ideal hydrodynamics for the system. 
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More term should be included in ([5]) once viscosity is The momentum rapidity of the flow 
considered. 

The proper time coordinate t, space-time rapidity co- 
ordinate 77, r and cj) are often used to note each space-time 
point. Thus 



yflo 



2 U(i-U^ 



(11) 



X = r cos 0, 

y — r sin0, 

z = T sinhry, 

f = r cosh?/. 



Or inversely 



1, t + z 
— in . 

2 i - z 



(3) 



(4) 



The initial state of the Landau hydrodynamical model 
is a thin static disk of thickness A and diameter a. It is 
assumed to be the approximation of the overlap of two 
highly Lorentz contracted nuclei, and the two geometrical 
sizes are assumed to be related by 



A = a/7, 



(5) 



with 7 the Lorentz factor corresponding to the velocity 
of the colliding nuclei in the center-of-mass frame. 
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VS. 



NN 



2m 



N 



(6) 



where Smn is the invariant collision energy per colliding 
nucleons of mass m n ■ 

The Landau hydrodynamical evolution of the system 
is divided in two parts: the initial longitudinal expan- 
sion (during which transverse velocities and displace- 
ments are neglected), and the subsequent "conic flight", 
where transverse velocities appear and hadronic particles 
freeze out of the fluid. 

Assuming a simple dynamical EOS of the form 



P = e/3, 



(7) 



an approximate solution of the 1+1 dimensional prob- 
lem of the longitudinal expansion phase is given by Lan- 
dau . This solution is summarized (using the notation 
of 0) a.s follow. The energy-density field is 



eo cxp 



{y+ + y-- Vv+y-) 



where y± are logarithmic light-cone coordinates 

, t±z 



(8) 



(9) 



and eo is the initial energy-density of the disk. The flow 
four- velocity are expressed as 



(cosh?7, 0, 0, sinhr/). 



(10) 



coincides with the space-time rapidity coordinate 77 de- 
flned in Eq.Q, similarly to the case of the Bjorken model. 

Now is the step of "conic flight", with hadronic par- 
ticles freezing out of the fluid. The four-momentum of 
each particle can be specified as 



= [jTiT cosh y,pT cos (f)p,pT sin (f)p,mT s'mhy) (12) 

in terms of its transverse momentum pT, rapidity y and 
azimuth angle 0p, where rriT — \/ m? + p^ is the trans- 
verse mass, with the mass of the particle, tti. 

The distribution of emitted hadrons in the momentum 
space is calculated with the Cooper-Frye formula Q 



dycPpT 



f{x,p)pt'da^, 



(13) 



where f{x,p) is the phase-space distribution of particles 
at freeze-out. The phase-space distribution has the ther- 
mal form 



fix,p) = f 



~ ( p'^Ufj,{x) 



T{x) 
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1 



(27r)3 



-P(^)+C' 

(14) 

with ^=1 for fermions, -1 for bosons and for classi- 
cal particles obeying Boltzmann statistics, g is a degen- 
eracy factor that is different for each particles species. 
The temperature at a given space-time point T{x) is ob- 
tained from the energy-density field in eq.([5]) with a cer- 
tain EOS. 

In previous work a hadronic gas EOS has been 
used, where the massive particles in the hadronic phase is 
assumed to obey the relativistic Boltzmann distribution. 
Then the pressure of the medium is given by 



P{T)^Y.'-^,mlT^K,{m^/T), 



(15) 



where rria and Qa are the mass and degeneracy factor of 
particle type a, and K2 is a modified Bessel function of 
the second kind. The energy-density and temperature is 
related via 



e(T)=3^P„ (1 

a 

where we used the notation 



mg Kia 
ST K20 



K,{mg/T). 



(16) 



(17) 



In Fig. [1] this relation is plotted with dashed line and 
note as "Pion", since the lightest mesons contribute the 
most to the pressure and energy density of hadronic gas. 
The same relation calculated with quark and gluon de- 
gree of freedom has been plotted with solid line and 
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Figure 1: (Color Online) Temperature dependence of en- 
ergy density. Dashed line; calculated according to Pion gas; 
Solid line: calculated with QGP freedom; Dotted-dashed line; 
LQCD result 0. 



noted as "QGP". We also plot Lattice result 01 as dotted 
dashed line and note as "LQCD". One can see Lattice 
EOS coincides with QGP one at high energy density and 
coincides with Pion one at low energy density. It is clear 
that no pions at very high energy density and no QGP 
at very low temperature, and LQCD result is a more re- 
alistic choice. In the following we repeat the previous 
calculation with Pion EOS and also take LQCD EOS as 
a comparison. 

We will employ Landau fluid dynamics in to RHIC and 
LHC heavy ion collisions, in the most simple case, cen- 
tral collisions. Therefore, the Landau solution of energy 
density in eq.([8|) is plotted as a function of rj and r in 
left panel of Fig. [21 where upper panel for the central 
Pb+Pb collision at ^s^vtv = 2.76 TeV and lower panel 
for the central Au+Au collision at ^/s^Jf = 200 GeV. 
The corresponding temperature is plotted in the right 
panel, based on two choices of EOS. 

To calculate hadrons' spectrum with Cooper-Frye for- 
mula in eq. (|13[) . we treat the hypersurface element four- 
vector 

(icr^i = {dxdydz, dydzdt, dxdzdt, dxdydt) (18) 

at each space-time point of FO surface, in terms of r, 77, r 
and (p according to eq. ([3]) for the convenience of calcula- 
tion. 

Insert u'^Uf^ = 1 into eq. (fT3|l . and make use of eq. ([T0| . 
one can get 



p^M^i — rriT cosh(y — rj) 



and 



'dau 



-TFodr], 



(19) 



(20) 



where the integral over r and </> has been done. 

In Refs. 0, |j| an estimate of the time of freeze-out is 



We take this as a choice of FO condition. Currently the 
popular FO condition is at a fixed temperature around 
Tc ^ 160 MeV. Here we take a constant energy den- 
sity epo a-s another choice of FO condition, which cor- 
responds to a fixed temperature via EOS. Since FO is 
a change of the treatment on collision systems, from a 
collective mode to the collection of many single parti- 
cles, this change should happen after the appearance of 
hadrons. Therefore, a realistic FO condition, should not 
occur at too high energy density where the matter is still 
in QGP phase. According to Fig. 1, three values 0.08, 
0.16 and 0.24 GeV/fm'^ were employed for epo, because 
Pion EOS and LQCD EOS coincide at e - 0.08 GeV/fm^ 
and 0.24 GeV/fm^ roughly corresponds to Tc ^ 160Mey 
according to LQCD EOS. 
With the FO condition 



(22) 



and the solution of Landau fluid in eq. , we can get the 
freeze-out time at each given space point as the following 
form, 



Tpo = 2a. 



(21) 



TPO = A exp ^ '- , 23 

6 

where the constant B = |ln(eo/epo), eo the initial 
energy-density of the disk, A the initial thickness of the 
disk and j] the space-time rapidity of the given space 
point. Then the integral along freez-out hypersurface 
can be done easily with Eqs. (fT9l [20]) . 



III. RESULTS AND DISCUSSION 

We first repeated the previous work [Sj], calculated the 
rapidity spectra of charged pions from central Au+Au 
collision at ^/sWn = 200 GeV with the same FO condi- 
tion Tpo = 2a and the same EOS. The result is shown 
as a dashed line in Fig. 3, a very nice coincidence with 
BRAHMS data where empty squares for tt^ and empty 
cycles for tt". However, this nice coincidence is just an 
accidence. We explain this in the following. 

In Fig. 3, the dashed line was decomposed into several 
dotted lines which are the contributions from different 77 
bin [-6,-4], [-4,-2], [4,6[, respectively. The energy den- 
sity and the temperature along 77-axis at rpo = 2a cor- 
respond to the second lines on the lower panel of Fig. 2. 
As one can see, most midrapidity pions come out from 
the fluid at 77 where the temperature is as high as 
300 MeV. This is not realistic because hadrons should 
not appear at such a high temperature, where the bary- 
onic chemical potentil is almost zero. In fact, this energy 
density is too high to employ Pion EOS, according to 
Fig. 1. 

Then we took LQCD EOS instead. This decreases the 
temperature at rpo = 2a quite a lot, c.f. the right panel 
of Fig. 2. In Fig. 4, the resulted rapidity distribution of 
charged pions and charged Kaons, shown as solid lines. 
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Figure 2: (Color Online) Left panehenergy density dependence of 77 at different proper time r where upper panel for central 
Pb+Pb collision at y^sjviv ~ 2.76 TeV and lower panel for the central Au+Au collision at -ysjvjv ~ 200 GeV. Right panel: the 
corresponding temperature as a function of 77 and r with two choices of EOS. 
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Figure 3: (Color Online) Rapidity distribution of charged pi- 
ous in central Au+Au collision at ^/snn ~ 200 GeV. Dashed 
line is calculated with Pion EOS and decomposed dotted lines 
for the contributions from r; bin [-6,-4], [-4,-2], 14,6], respec- 
tively. Data points are from Q. 




and Pion EOS results shown as dashed lines, are com- 
pared with BRAHMS data points 0]. While Pion EOS 
results can reproduce BRAHMS measured pion and Kaon 
rapidity distribution, this realistic LQCD EOS calcula- 
tion can not! 

The reason is clear. At rpo — 2a, hadrons of large 
rapidity mainly freeze out from the fluid at large- 77, where 
energy density is low and LQCD EOS coincides with Pion 
EOS. While at midrapidity, hadrons mainly freeze out 
from fluid at 77 ^ where the energy density is high. 
Then the temperature interpreted with LQCD EOS is 
much lower than that with Pion EOS, which makes less 
hadron emission at midrapidity. As a result, the rapidity 
shape with LQCD EOS becomes broader and deviates 
from the measured data. The correct choice of LQCD 
EOS itself can not reproduce the data, but the deviation 



Figure 4: (Color Online) Rapidity distribution of pion(upper 
panel) and Kaon(lower panel) in central Au+Au collision at 
y^sWN = 200 GeV, with Pion EOS(dashed line) and LQCD 
EOS(solid line). Data points are from 0]. 



can tell us even more. 

To understand the deviation better, we calculated the 
transverse mass spectra of charged pions and Kaons from 
central Au+Au collisions at y^s^Jf ~ 200 GeV, with the 
FO condition rpo = 2a and two options of EOS. The 
result with Pion EOS is shown in Fig. 5, where solid 
lines for pions and dashed lines for Kaons, at midrapid- 
ity y = (left panel) and forward rapidity y = 3.5 (right 
panel) and compared with BRAHMS data where full cy- 
cles for n~ and full squares for K~ . The calculation with 
LQCD EOS is plotted in Fig. 6, where the same notation 
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has been used. Again, similar results appear at forward 
rapidity y = 3.5 for both EOS, where particles freeze 
out from quite low energy density source and two EOS 
coincidence. 

However, at midrapidity y = 0, the inverse slope of the 
transverse spectra of both pions and Kaons from LQCD 
calculation is lower than experimental data, while the cal- 
culation with Pion EOS makes higher the inverse slope 
for pions, due to the extremely high temperature inter- 
preted with Pion EOS. 

The second line in lower-right panel of Fig. 2 tells the 
temperature interpreted with LQCD is about 200 MeV 
at 77 ^ 0, which is lower than Pion EOS value, but still 
higher than Tc ^ 160 MeV. Then Why the calculated in- 
verse slopes of both pions and Kaons are still lower than 
experimental data? Because there is no transverse expan- 
sion in Landau hydrodynamics. Hydrodynamical models 
with transverse expansion and freeze-out at 160 MeV can 
make higher inverse slope of transverse spectra and repro- 
duce the data, i.e. The deviation tell us, the collective 
transverse expansion of the collision systems is important 
to explain the measured transverse spectra. 

The midrapidity transverse momentum spectra of 
charged pions, Kaons and antiproton from central 
Pb+Pb cohisions at = 2.76 TeV, with the FO 

condition rpo = 2a and two options of EOS are shown in 
Fig. 7, where solid lines for Pion EOS and dashed lines 
for LQCD EOS, and ALICE data O shown as fuh cy- 
cles. Again we can see Pion EOS gives high inverse slopes 
than LQCD. The calculated pion slope with LQCD EOS 
seems to reproduce the data, again by accidence. Be- 
cause the FO condition rpo = 2a makes a much higher 
energy density and the temperature at freeze-out is in- 
terppreted as high as 400 MeV, c.f. the second line in the 
upper panel of Fig. 2. This compensates the missing of 
transverse expansion in Landau fluid. 

In Fig. [71 the deviation of transverse spectra between 
data and LQCD results, becomes bigger and bigger, as 
the mass ordering from pion, Kaon to antiproton. This is 
because, transverse flow drives heavier particles to move 
along radial direction at the same collective velocity as 
light particles. However, heavier particles gain more en- 
ergies and higher inverse slope of the transverse momen- 
tum spectra. In another word, transverse expansion of 
the collision systems becomes more and more important 
with the increase of hadron mass. 

As mentioned above, the FO condition rpo ~ 2a is 
not realistic, because the energy density and temperature 
vary a lot along the beam direction. In the following, we 
take a constant energy density epo ^-s a FO condition. 

In Fig. 8 are shown the transverse mass spectra of pion 
(upper panel) and Kaon (lower panel) in central Au+Au 
collision at ^/Tnn ~ 200 GeV, where left panel for y=0 
and right panel for y=3.5, with solid line, dashed line 
and dotted line for eFo=0.08, 0.16, 0.24 GeV/fm^, re- 



spectively. The higher epo can slightly increase the in- 
verse slope of the transverse spectra, while the slopes 
are still far from the BRAHMS data, due the missing of 
transverse flow in Landau fluid. 

In Fig. 9 are shown the midrapidity transverse mo- 
mentum spectra of charged pions, Kaons and antiproton 
from central Pb+Pb collisions at ^fs^N = 2.76 TeV, 
with the FO condition epo=0.08, 0.16, 0.24 GeV/fm^ 
shown as solid lines, dashed lines and dotted lines re- 
spectively. ALICE data [l^l are shown as full cycles. 
Again, higher epo makes slightly higher inverse slope. 
But compared with Fig. 8 for the case of central Au+Au 
collision at ^/snn = 200 GeV, the miss of transverse ex- 
pansion of Landau fluid in central Pb+Pb collisions at 
y/sNN ~ 2.76 TeV is even stronger, for the larger and 
hotter collision system. 



IV. CONCLUSIONS 

In summary, we demonstrated the production of 
most copious particles from central Au+Au collisions 
at y^SNN = 200 GeV and central Pb+Pb collisions at 
y/SNN =2.76 TeV in the framework of Landau hydrody- 
namical model. 

The measured rapidity distribution of pions and kaons 
from central Au+Au collisions at y^SNN = 200 GeV and 
the transverse momentum spectrum of pions from Cen- 
tral Pb+Pb collisions at ^/snn = 2.76 TeV can be re- 
produced with Landau model by accidence, due to the 
choice of FO condition rpo = 2a and the employment 
of Pion EOS. This FO condition makes most hadrons 
freeze out from a fluid at very high energy density. The 
employment of Pion EOS causes that the temperature at 
which hadrons are emitted is interppreted to be as high 
as 400 MeV. The emission of hadrons at high density and 
high temperature compensates the missing of transverse 
expansion in Landau hydrodynamical model. 

With a more realistic choice of FO condition, fixed 
at a relatively low energy density, and the employment 
of LQCD EOS for particle production, we can see that 
transverse expansion becomes important to explain bulk 
hadron production in high energy nuclear collisions. Es- 
pecially, with the increase of collision energy from RHIC 
to LHC, and with the increase of the mass of produced 
hadrons, transverse expansion becomes more and more 
important. 
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Au+Au \[sMfj=200GeV 0-5% x = 2a, Pion EoS 
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Figure 5: (Color Online) Transverse mass spectra of charged pions and Kaons from central Au+Au collisions at ^sjvjv = 
200 GeV at midrapidity y = (left panel) and forward rapidity y = 3.5 (right panel). The calculation was done with Pion 
EOS and FO condition rpo = 2a, where solid lines for pions and dashed lines for Kaons. Cycles for n~ and squares for K~ of 
BRAHMS data Q. 



Au+Au \[sM|;j=200GeV 0-5% x = 2a, LQCD EoS 




m-r-ni[GeV/c^] mT-nj[GeV/c^] 



Figure 6: (Color Online) Same as Fig. 5 except the calculation was done with LQCD EOS. 



Pb+Pt\'sMM=2.76TeV 0-5%x,„=2a 




Figure 7: (Color Online) Transverse momentum spectra of pion ,Kaon and anti-proton from central Pb+Pb collisions at 
^/snn ~ 2.76 TeV. FO condition rpo — 2a, solid lines for Pion EOS and dashed lines for LQCD EOS. Data points from ALICE 
collaboration [To| . 
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Au+Au\[S;;jN=200GeV 0-5% LQCD EoS 




Figure 8: (Color Online) The transverse mass spectra of pion (upper panel) and Kaon (lower panel) in central Au+Au collision 
at y'SNN ~ 200 GeV, where left panel for y=0 and right panel for y=3.5. Solid lines, dashed lines and dotted lines for FO 
condition epo=0.08, 0.16, 0.24 GeV/fm^, respectively. Data points are from the BRAHMS collaboration 
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Figure 9: (Color Online) Transverse momentum spectra of pion ,Kaon and anti-proton from central Pb+Pb collisions at 
y/SNN ~ 2.76 TeV. Solid lines, dashed lines and dotted lines for FO condition eFo=0.08, 0.16, 0.24 GeV/fm^, respectively. 
Data points from ALICE collaboration [lo| . 
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